arXiv:1502.04976vl [astro-ph.GA] 17 Feb 2015 


Tully-Fisher Scalings and Boundary Conditions for 

Wave Dark Matter 


Andrew S. Goetz 
February 18, 2015 


Abstract 

We investigate a theory of dark matter called wave dark matter, also known 
as scalar field dark matter (SFDM) and boson star dark matter or Bose-Einstein 
condensate (BEC) dark matter (also see axion dark matter), and its relation to 
the Tully-Fisher relation. We exhibit two boundary conditions that give rise to 
Tully-Fisher-like relations for spherically symmetric static wave dark matter halos: 
(BCl) Fixing a length scale at the outer edge of wave dark matter halos gives 
rise to a Tully-Fisher-like relation of the form = constant. (BC2) Fixing 

the density of dark matter at the outer edge of wave dark matter halos gives rise 
to a Tully-Fisher-like relation of the form = const. These results extend 

the work of a previous paper . 


1 Introduction 


The baryonic Tully-Fisher relation (BTFR) is an empirical scaling relation for disk 
galaxies involving the baryonic (i.e., non-dark) mass M}, and the rotational velocity v 
of stars. It is observed that 

Mb 


— = const. 

yX 


( 1 ) 


where the exponent x (the slope of the relation in log-log space) is usually found to 
lie between 3 and 4. There is an extensive literatnre which nses ristronomical data to 
calibrate the relation and estimate x. See, for example, [I§|Il|2g|l6l@i|I||27l|l7l. 

In the relation ([^, the rotational velocity v is determined by the total gravitating 
mass (baryonic and dark) while acconnts for the baryonic matter only. Thns, the 
BTFR provides a link between the baryonic and dark matter and may have something 
to say abont the viability of varions theories of dark matter. 

In this paper we discuss the BTFR in the context of a theory of wave dark matter. 
Wave dark matter 20, 21 has been investigated under other names such 


as scalar held dark matter (SFDM) 


I0|[^ [l4l[26 


22 


and boson star dark matter or 
Bose-Einstein condensate (BEC) dark matter (^|^[l2, The difference 


in names comes from a difference in motivations, but the underlying equation is the 


Klein-Gordon wave eqnation (2b) for a scalar held. 
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In spherical symmetry, there are simple “static states” of wave dark matter (see 
hgs. to at the end of the paper). Sin proposed in 25 that each galactic dark 
matter halo corresponds (in a hrst approximation) to one of these static states. If 


we take this snggestion as a starting point, a natnral concomitant idea is that the 
baryonic Tully-Fisher relation might arise because a corresponding “Tully-Fisher-like” 
relation holds for the wave dark matter halos. This Tully-Fisher-like relation would 
hold because of the nature of wave dark matter. The fundamental question we ask 
in this paper is whether the equations of wave dark matter have the potential to give 
rise to such a Tully-Fisher-like relation (see Question on page|^. We answer this 
question in the affirmative and exhibit two boundary conditions, which could be physical 
conditions imposed at the edge of each wave dark matter halo, that give rise to just 
such a Tully-Fisher-like relation. The boundary conditions are as follows: 


BCl: Fixing a length scale at the outer edge of halos implies a Tully-Fisher-like 
relation with slope 4. This is a more general version of a result from a previous 
paper [^. 

BC2: Fixing the dark matter density at the outer edge of halos implies a Tully-Fisher- 
like relation with slope 3.4. 

The remainder of this paper runs as follows: In section we introduce the spherically 
symmetric static states and their scaling properties. In section we exhibit BCl 
and BC2 numerically, and in section we derive BCl and BC2 theoretically. In 
section we describe other conditions one might impose which do not lead to Tully- 
Fisher-like relations. Finally in section we briefly discuss the connection between 
a Tully-Fisher-like relation for wave dark matter halos and the baryonic Tully-Fisher 
relation. 


2 The Spherically Symmetric Static States 


The equations underlying the theory of wave dark matter are the Einstein-Klein-Gordon 
(EKG) equations: 


G 


□»/ 


Stt 


/ df ®df + df ®df 
T2 


TV- 



(2a) 

(2b) 


(Throughout, we use geometrized units in which the universal gravitation constant and 
the speed of light are unity: G = c = 1.) Here / is a complex scalar field representing 
the dark matter and T is a fundamental constant of nature. The relationship between 
T and the wave dark matter particle mass m is 

m = — = 2.09 X 10-23 eV —^ . 3 

c \l\y V 
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In this paper for definiteness we use T = 101y“^, which is consistent with constraints 
from other studies . 

Because we are interested in studying galactic dark matter halos, most of which are 
thought to be approximately spherically symmetric, we study the EKG equations in 
spherical symmetry. The spacetime metric we will use is 

g = dt^ + r)~^ + r‘^{d6‘^ + siv? (f) dcf)^) (4) 

where $(t, r) = 1 — the non-relativistic limit, the functions V{t,r) and 

M(t, r) have natural interpretations as the gravitational potential at time t and radius 
r and the total mass contained inside a ball of radius r at time t. The simplest solutions 
to eq. (|^ in the metric Q are the spherically symmetric static states, which are of the 
form 

/(f,r)=F(r)G"b (5) 

Substituting the ansatz (|^ into eq. ([^ and using eq. Q, we obtain (see HD) the 
relatively simple system of ODEs 

Mr = dvrr" • ^ [(T^ + (6a) 

^ F^ - ^F^] (6b) 

Frr + ‘^Fr + (T^ - F (6c) 

which can be integrated on a computer. Here the functions M{r),V{r),F{r) are 
functions of r only, and the only dependence on t appears in eq. ([^. 

In the non-relativistic limit the solutions of ([^ are indistinguishable from the 
solutions of the even simpler system of ODEs 

Mr = Airr^ ■ 2F^ (7a) 

M 

(7b) 

^ (^Frr + ‘^-F^j ={T-UJ + TV)F (7c) 


This system is the Poisson-Schrddinger (PS) system written in spherical symmetry. Here 
we have a special case of the well-known fact that the PS system is the non-relativistic 
limit of the EKG system—see . From here on in this paper we work exclusively with 
the simpler PS system ([^ because we are interested in wave dark matter in galaxies 
which are, considered as an aggregate, non-relativistic systems. 

Finite-mass solutions to eq. ([^ come in the form of ground states and excited states. 
It was proposed in 25 that each galactic dark matter halo corresponds (in a hrst 
approximation) to one of these states. In hgs. [^to|^we display some representative 
graphs of solutions so that the reader can get a sense of what they look like, and refer 
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the reader to a previous paper for more examples. We assign each solution a non¬ 
negative integer n which counts the number of zeros (nodes) of F. A ground state has 
n = 0; excited states have n ^ 1. Any solution can be given as a quadruple (cn; M, V, F) 
consisting of the constant u followed by the three functions M{r), V"(r), and F{r). It is 
easily verihed that if (cu; M, V, F) is a solution, then so is {u + TV] M, V + V,F). We 
use this freedom to shift every potential function V{r) so that it approaches zero at 
infinity. 

For each static state there is a special value of r which we refer to as Rbm and 
which we think of as the edge of the dark matter halo. Mathematically, Rbm is 
the particular value of r at which the function F{r) switches from oscillatory to 
exponentially decreasing behavior. Looking at eq. (7c), we see that Rbm satisfies the 
equation T — a; -I- Tl/(i?DM) = 0. Thus we make the following definition: 


Definition 1. Given any spherically symmetric static state we define 7 ?dm to be the 
radius at which the function F switches from oscillatory to exponential behavior. It 
satisfies the equation 

T - a; + TG(i?DM) = 0. (8) 

We regard r = Rbm as the edge of the dark matter halo because almost all the dark 
matter is contained in the region r ^ Rbm- See figs. |^to[^ 

It is well-known that solutions to the PS system ([^ have certain scaling properties. 
Specifically, for a,/3 > 0, we can scale any given solution as follows: 


r = 


(9a) 

M = 

a^-^M 

(9b) 

V = 

a^/3-^V 

(9c) 

F = 

a^F 

(9d) 

T = 

(3‘^T 

(9e) 

(t-w). 

Q:^(T — Oj) 

(9f) 

se scalings 

can be derived immediately from eqs. 

(7a) 


to (7c). Note that (3 effectively controls the fundamental constant of nature T and 
a can be used to control the “peakedness” of the solution. States which are more 
compact have higher mass; states which are more spread out have lower mass. In fig. 
we show the effect of scaling a ground state using values of a near 1. It follows from 
these scalings that having hxed T, there is a one-parameter family of static states of 
order n. Once we fix a further characteristic of the solution, a static state is uniquely 
determined for each n, giving us a sequence of static states. We call this “hxing a 
scaling” or “imposing a scaling condition”. 
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Figure 1: Five scalings of the n = 1 excited state from fig. (shown in blue) using 
Oi = 1.1,1.2,1.3,1.4,1.5 in the scaling formulas ([^. Here we only show the dark matter 
prohle function F(r). 


3 Tully-Fisher Scalings and Boundary Conditions: 
Numerical Evidence 


The basic question we are asking in this paper is: 


Question 1. Are there are any scaling conditions for which the sequence of static 
states n = 0,1,2,... obeys a Tully-Fisher-like relation 


M 

— = constant 

vx 

for some 3 < a; ^ 4? 


( 10 ) 


We need to make this question more precise, as it is unclear above what M and v are, 
exactly, given a particular static state. Astrophysicists attempting to calibrate the 
baryonic Tully-Fisher relation are faced with a similar problem. In that case, M is an 
estimate for the total baryonic mass. For v, some take it to be the maximum observed 
velocity; some take it to be an average velocity computed in some well-dehned way; 
some take it to be the velocity observed at a particular (arbitrarily chosen) radius. It is 
well-known that the rotation curves of many spiral galaxies are somewhat flat, and in 
this case any of these three choices will produce similar values for v. 

We have chosen to take M = M(i?DM) and v = u(i?DM)- This latter quantity is 
the circular velocity which would be observed for objects orbiting at r = Rbm- This is 
computed in the standard Newtonian way: 


u(r)^ M(r) 

lyt rjr^‘2^ 



( 11 ) 
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Perhaps a more natural choice for M would be M^o = limr^oo M{r), i.e. the total mass 
of the static state. But M^Rdm) Moo since Rbm is by dehnition the “outer edge” 
of the halo (see hgs. |^to|^, so it does not really matter which we choose. Similarly, 
another natural choice for v would be Umax, the maximum value of the function v{r). 
But again, by examining figs. to one can see that v{Rdm) ^ U ma x. Thus, the 
following more precise version of Question seems reasonable: 

Question 2. Are there are any scaling conditions for which the sequence of static 
states u = 0,1,2,... obeys a Tully-Fisher-like relation 

M^Rbm) ^ ('X2) 

v{Rdm)^ ^ ^ 

for some 3 < a; < 4? 


Now, in one sense the answer to Questionis a trivial “yes”, for we can just pick 
any x we like and let eq. (12) itself be the scaling condition. However, this is obviously 
cheating; we are interested in scaling conditions which might correspond to something 
physical. In fact, we have found two conditions that we consider to be more physical. 
These were given in section along with the Tully-Fisher-like relations thay imply. We 
restate them here more precisely: 

BCl: Fixing a length scale at Rbm implies a Tully-Fisher-like relation with slope 4. 


BC2: Fixing |T(i?DM)| implies a Tully-Fisher-like relation with slope 3.4. 

With regards to BCl, we should explain what “fixing a length scale” means. We 
postpone a full explanation to section and limit ourselves here to an explanation 
which involves what we term the “halflength” of a state, defined as follows: 

Definition 2. The F{r) function for any static state has an exponentially decreasing 
“tail” lying to the right of i?DM- There is some R > Rbm such that F{R) = ^F{Rbm)] 
we refer to R — Rbm as the halflength of the state. 


Requiring the halflength of all states to be the same is one way to hx a length scale at 
Rbm- We refer the reader to section for full details. 

With regards to BC2 note that, by eq. (7a), |F| determines the density of the wave 
dark matter. This explains why in section [T we referred to “hxing the dark matter 
density” in BC2. 

It remains in this section to show numerically that BCl and BC2 really do give 
Tully-Fisher-like relations with the stated slopes. We used a computer to derive 
particular static state solutions for n = 0,1,..., 800 and then scaled these solutions in 
accordance with eq. 0 so that, first, they all had the same halflength, and second, so 
they all had the same value for |T(i7DM)|- (The precise values for the halflength and 
|T(Rdm)| were chosen so that the resulting static states would be of reasonable sizes 
to be representing galactic halos.) The results are shown in hg. in which we have 
plotted M{Rbm) vs. ^(Rdm) m log-log space. For large n, the static states form a line 
with the stated slope. (For small n, the states deviate slightly from the correct slope. 
This will be explained in section]^) 
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Fixed Halflength 


Fixed |F(/2dm)| 




Figure 2: Numerical evidence for BCl and BC2 (see page|^. At left, the static states 
n = 0,1,... ,800 with T = 101y“^ and a fixed exponential decay halflength of 884 ly. 
The dashed line is the has slope 4. At right, the static states n = 0,1,... ,800 with 
T = 101y“^ and a fixed value for |F(i?DM)| of 2.473 x 10“®. The dashed line has slope 

3.4. 

4 Tully-Fisher Scalings and Boundary Conditions: 
Theoretical Derivations 


In this section we wish to explain theoretically why the statements made in BCl and 
BC2 are true. We have already exhibited their truth numerically in fig. 

Both statements involve boundary conditions which fix some aspect of the solutions 
to the PS system ([^ around the point r = Rbm- We can approximate the solutions 
around i?DM by using a linear approximation for V{r): 


V{r) ^ + P'(i?DM)(r - Rdm) 

M (Rdm) 


— y (Rdm) + 


C)2 

-'^DM 


(r — R 


DMJ 


for r ss Rbm- Substituting this approximation for V into eq. (7c) and using eq. 
And that for r % Rbm, F{r) approximately solves the differential equation 

^ , 2^ 2T^M{Rbm),^ „ 

Frr T Fj- 2 y RbMjF. 

r Rhu 


(13) 


we 
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Multiplying through by r, we can write this differential equation as 

2T‘^M{Rdu) r> \r Tr\ (^a\ 

{rF)„ = -- [r - R^u){rF). (14) 

-^DM 

Compare this to the Airy differential equation y" = xy whose general solution is a linear 
combination of the two standard functions Ai(a:) and Bi(x). Both Ai(x) and Bi(x) 
oscillate for negative x and have exponential behavior for positive x. The function Ai(x) 
is shown in fig. It exponentially decreases to zero for x > 0, whereas the function 
Bi(x) (not shown) exponentially increases for x > 0. 



Figure 3: In black, the standard solution Ai(x) to the Airy differential equation y" = xy. 
In brown and dashed, the approximation (21) to Ai(x) for x < 0. 


Looking at eq. (14), we see that for r % rF approximately solves an equation 

of the form y" = a{x — xo)y. The general solution to this differential equation is 
A ■ Ai(a^/^(x — xo)) + B ■ Bi(a^/^(x — xq)). Since F exponentially decreases for r > Rbm, 
we must have, therefore. 


F(r) oc - Ai 


2T^MiRj,u)V^" 


-^DM 


-Rdm) 


(15) 


for r ss i?DM- If we choose the proportionality constant in eq. (15) so that the 
approximation has the correct value at r = Rbm, then we get 




2TUf(flDM) 

-^DM 


1/3 


(r — Rbm) 


(16) 


for r ss i? 


DM- 


Equation (15) demonstrates mathematically the truth of BCl. For up to a vertical 


scaling, the shape of F{r) around r = Rbm is given by the expression on the right side 


of (|15|). Thus, fixing a length scale around r = i?DM evidently amounts to fixing the 

2T^M{R^m) 


value of the constant 


-^DM 
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(17) 



















This argument explains why fixing a length scale is basically equivalent to fixing a value 
for the halflength (as we did in section or to any of a nnmber of other conditions 
which fix some length aronnd -Rdm (for another example, see onr previons paper 0 ), 
All such conditions are fixing a particular horizontal scaling of the expression in eq. (16), 
Now, the rotation curve calculation ([IT|) shows that 




jd2 

■^DM 


(18) 


since T is a constant (in this paper, T = lOly ^), we see that fixing a length scale is 


equivalent to fixing the constant in eq. (18). This explains why we get the Tully-Fisher- 


like relation shown in fig. [^with slope 4. The approximation (15) gets better for larger 
n (see figs. |^to[^, which explains why in fig. [^the lower order static states deviate 
from the slope 4 line. 

To show the truth of BC2, we begin with the approximation (16). This expression 
contains the four numbers T, Rdm, M{Rdm), and F{Rdm)- These cannot all be 
independently chosen, for the static states effectively live in a 3-parameter space. 
(Recall the 2-parameter family of scalings (|^. A third parameter hiding there is the 
discrete parameter n, the order of the static state.) Thns, we shonld be able to find 
some relation between these four numbers. In fact, we will argue that 

|F(Rdm)| ^ Ai(o) . for large n (19) 

where 

2^17/12 Ai(o) ^ 0.133 (20) 

This explains BC2, for then we see that fixing \F{Rj)m) \ amounts to fixing M(i?DM)^-RDM' 
By eq. (jm, Rdm = M{R-[)m)v{Ry)m)~‘^, so fixing M(Rdm)^Rdm is equivalent to fixing 
fof(RDM)^°n(RDM)^^ = {M{Ry)m)/v{Rt)m)^'^)~^^, i-e. equivalent to fixing 

M (Rdm) 

^(Rdm)^'^ 


Our goal now is to derive eq. (19) beginning with the approximation (16). We will 


utilize a well-known approximation for the Airy function Ai(x) for negative inputs: 

( 21 ) 


sin -t-f) 

Ai(—a:) % - —— for x > 0. 


TTX 


1/4 


This approximation can be seen in fig. It is extremely good for x > 1. We will also 
use the following definite integral, which is easily obtained: 


r 


2 (2 3/2 , vrl 

^ 3 "^ ^ ildx = L^/^ + 0{l). 


X 


1/2 


( 22 ) 


The notation 0(1) is mathematical shorthand for a bounded function of L. 


9 













Using eqs. (7a) and (16), we have 

"’^DM 

M(i?DM) = I Svrr^F^dr 


I 

I 




Svr 


Ai(0)2 


Ai(0)2 


Ai 


2T2M(i?DM) 

-^DM 

2T2M(i?DM) 

d2 

-'^DM 


1/3 


n 2 


(-Rdm — r ) 


dr 


1/3 


H 2 


dr. 


We will temporarily abbreviate -Rdm = R, M{Rdm) = M, F(i?DM) = -U, Ai(0) = A. 

( 2 \ 1/3 

Making the subsitution x = ( and using the approximation (21), we obtain 


F2/?2 


M ss Svr- 


A2 


2T2]\^\-V3 p(2T2MR)i/3 


i?2 


28/3 

~A^ 


Y-2/3jp2j^8/3jy-l/3 


I 


Ai(—a:)^ dx 

[2T2Mii]V3 ^.^2 (2^3/2 + 


X 


1/2 


dx. 


Using the dehnite integral (22), we obtain 

28/3 


M 


A2 

217/6 


A2 


j-2/3p2j^8/3^-l/3 [(2T2Md2)^/® + 0(1)] 
'Y'-V3^2^17/6j^^-l/6 


(We dropped the 0(1) term because the other term dominates for large n.) Rearranging, 
we get 

F2 ^ 2-^'^/^A^T^/3m'^/^R-^'^/\ 


Taking square roots gives us eq. (19), which completes the argument. 

One might wonder how good the approximation (19) is for small n. In hg. 
have graphed the ratio 

|F(Rdm)| 


we 


2-l’/12 Ai(0) ■ 


(23) 


as a function of n. For the ground state (n = 0), the ratio is approximately 1.15, which 
is reasonably close to 1, and for large n the ratio is extremely close to 1. 

Combining eqs. (16) and (19) gives the following approximation for the F{r) function 
of a static state of order n. (By convention F{0) is positive so the factor (—1)” is 
necessary.) 

F(r) ^ (-1)"2-'^/'2^2 /i2^^^^^)7/i 2^-V12 , 1 f f2T^M{RuM)^ 


-^DM 


(r — Rdm) 


This approximation is shown in hgs. to 


(24) 
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5 Other Scaling Conditions 

The scaling conditions given in BCl and BC2 are boundary conditions imposed at the 
edge of dark halos. The hrst imposes a horizontal scale on F{r) and the second imposes 
a vertical scale. One might wonder whether there are any other scaling conditions which 
lead to Tully-Fisher-like relations. For example, what happens when we try the same 
thing at the center of the halo? To impose a horizontal scale on F{r) at the origin, 
we can fix the location of the hrst node (for n ^ 1). To impose a vertical scaling, we 
can hx |F(0)|. In hg. |^we show the results of imposing these two scaling conditions. 
Neither leads to a Tully-Fisher-like relation. In fact, besides BCl and BC2 no other 
scaling condition we have tried leads to a Tully-Fisher-like relation. 


6 Conclusion 

We have exhibited two boundary conditions, BCl and BC2, which we restate here 
once more in their more colloquial form from section 

BCl: Fixing a length scale at the outer edge of spherically symmetric wave dark 
matter halos implies a Tully-Fisher-like relation with slope 4. 

BC2: Fixing the dark matter density at the outer edge of spherically symmetric wave 
dark matter halos implies a Tully-Fisher-like relation with slope 3.4. 

We gave both numerical evidence and theoretical arguments illustrating the truth 
of BCl and BC2. Our main goal in this paper has been to demonstrate that a 
Tully-Fisher-like relation can arise naturally out of the equations of wave dark matter. 

To conclude we comment briefly on the connection between a Tully-Fisher-like 
relation for wave dark matter halos and the observed baryonic Tully-Fisher relation. 
Dehne the baryon fraction of a galaxy to be the fraction of the galaxy’s mass which 
is baryonic. If the baryon fraction were the same for every galaxy, then the baryonic 
Tully-Fisher relation would have the same slope as a Tully-Fisher-like relation for dark 
halos. However, it is observed that the baryon fraction of galaxies decreases with size. 
Small galaxies contain mostly dark matter, and only very large galaxies contain baryons 
in sufficient quantities to approach the cosmic baryon fraction, which is around 1/6. 
This is known as the missing baryon problem (see [I^). 

In M vs. V log-log space, when switching from plotting dark mass to baryonic mass, 
all galaxies shift vertically downward. The missing baryon problem implies that large 
galaxies shift less than small galaxies. Therefore, if the baryonic Tully-Fisher relation 
is related to a Tully-Fisher-like relation for dark halos, the former relation should have 
a steeper slope than the latter. 
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0 100 200 300 400 500 600 700 800 


Figure 4: The ratio given in eq. (23) plotted as a function of n (u = 0,1,..., 800). This 
graph shows that the approximation (19) is very good for large n and not too bad for 
small n. 


Fixed |T(0)| 


Fixed Location of First Node 



10 ^ 10 ^ 10 ^ 
v{Rbm) (km/s) 


0 


Q 

0 ^ 



^^(7?DM) (km/s) 


Figure 5: At left, the static states u = 1,..., 800 with T = 10 ly“^ and the first node 
fixed at 1500 ly. At right, the static states n = 0,1,... ,800 with T = 101y“^ and 
|F(0)| = 8 X 10“^. The dashed lines have slope 4. Neither scaling condition gives a 
Tully-Fisher-like relation. 
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Figure 6; A ground state {n = 0) with T = 101y“^. Units are geometrized and given 
in light-years. The first three functions graphed are the solutions M(r), V{r), F{r) to 
eq. 0 - They are, respectively, the total mass prohle, the Newtonian potential, and 
the dark matter prohle. The fourth function is a rotation curve v{r) calculated in 
accordance with eq. The vertical gray line marks the location of Rbm- The dashed 
curves show the approximations described in section]^ in eqs. (13) and (24). 


13 


















14 









8 

6 

4 

2 

0 

0 

,5 

1 

5 

2 

2 

1 

0 

2 

1 

0 


• 10 -^ 




• 10-7 




15 














Figure 9: An excited state {n = 100). See the caption to fig. [^for a fuller description. 
The top and bottom of the F{r) function have been cropped out to show more detail. 
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